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Abstract 

The melting of ultrathin lubricant film by friction between atomically flat surfaces 
is studied. The additive noises of the elastic shear stress and strain, and the tem- 
perature are introduced for building a phase diagram with the domains of sliding, 
stick-slip, and dry friction. It is shown that increase of the strain noise intensity 
causes the lubricant film melting even at low temperatures of the friction surfaces. 
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1 Introduction 



The study of the noise influence on the friction process has an evident fundamental and 
practical importance because in some experimental situations the fluctuations can change 
the frictional behavior critically, for example, providing the conditions for low friction 
[1] [3]. In particular, the thermal noise, acting in any experiments, can convert the 
ultrathin lubricant film from stable solidlike phase state to the liquidlike one and, thus, 
transform the dry friction into the sliding or the stick-slip (the interrupted) modes. There- 
fore, in recent years the considerable study has been given to the influence of disorder and 
random impurities in the interface on the static and the dynamic frictional phenomena 
[4] - - [6]. These investigations show that a periodic surfaces are characterized by smaller 
friction coefficient during sliding than nonregular ones. Besides, the stick-slip dynamics, 
inherent in solid friction, attracts an increased attention on the atomic [7] [9] and 
the macroscopic [10,11] levels as well as for granular mediums [12] [14]. In order to 
achieve the better understanding of the above phenomena, here an analytic approach is 
put forward, which describes the transitions between friction modes due to variation of 
fluctuations of elastic and thermal fields. 

In the previous work [15] on the basis of rheological description of viscoelastic medium 
the system of kinetic equations has been obtained, which define the mutually coordinated 
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evolution of the elastic shear components of the stress a and the strain s, and the temper- 
ature T in ultrathin lubricant film during friction between atomically flat mica surfaces. 
Let us write these equations using the measure units 



V t t J G \t t J \ T] ) 

for variables a, e, T, respectively, where p is the mass density, c v is the specific heat 
capacity, T c is the critical temperature, rjo = r](T = 2T C ) is the typical value of shear 
viscosity rj, tt = pl 2 c v / k is the time of heat conductivity, / is the scale of heat conductivity, 
k is the heat conductivity constant, r £ is the relaxation time of matter strain, G = r] /T £ : 



r a a = -a + ge, (2) 
T £ e = —e + (T — l)a, (3) 
r T f = (T e - T) - as + a 2 . (4) 

Here the stress relaxation time r CT , the temperature T e of atomically flat mica friction 
surfaces, and the constant g = G/Gq are introduced, where G is the lubricant shear 
modulus. It can be seen [15,16] that Eqs. (2) and (3) are the Maxwell-type and the Kelvin- 
Voigt equations for viscoelastic matter, correspondingly. The latter takes into account 
the dependence of the shear viscosity on the dimensionless temperature rj = i] /{T — 1). 
Equation (4) represents the heat conductivity expression, which describes the heat transfer 
from the friction surfaces to the layer of lubricant, the effect of the dissipative heating of 
a viscous liquid flowing under the action of the stress, and the reversible mechanic-and- 
caloric effect in linear approximation. These equations coincide with the synergetic Lorenz 
system formally [17,18], where the elastic shear stress acts as the order parameter, the 
conjugate field is reduced to the elastic shear strain, and the temperature is the control 
parameter. As is known this system can be used for description of the thermodynamic 
phase and the kinetic transitions. 

In Ref. [15] a melting of ultrathin lubricant film by friction between atomically flat mica 
surfaces has been represented as a result of action of spontaneously appearing elastic field 
of stress shear component caused by the heating of friction surfaces above the critical 
value T c = 1 + g -1 . Thus, according to such approach the studied solid-liquid transition 
of lubricant film occurs due to both thermodynamic and shear melting. The initial reason 
for this self-organization process is the positive feedback of T and a on e [see Eq. (3)] 
conditioned by the temperature dependence of the shear viscosity leading to its divergence. 
On the other hand, the negative feedback of a and e on T in Eq. (4) plays an important 
role since it ensures the system stability. 

According to this approach the lubricant represents a strongly viscous liquid that can 
behave itself similar to the solid — has a high effective viscosity and still exhibits a yield 
stress [7,16]. Its liquidlike state corresponds to the elastic shear stress which relax to the 
zero value during larger times r CT than that of the solidlike state. Moreover, at a = 
Eq. (3) reduces to the Debye law describing the rapid relaxation of the elastic shear strain 
during the microscopic time r £ ~ 1CT 12 s. At that the heat conductivity equation (4) 
takes the form of simplest expression for temperature relaxation that does not contain 
the terms representing the dissipative heating and the mechanic-and-caloric effect of a 
viscous liquid. Also, we assumed that the film becomes more liquidlike and the friction 



decreases with the temperature growth due to decreasing activation energy barrier to 
molecular hops. 



In accordance with Ref. [11] in the absence of shear deformations the temperature mean- 
square displacement is defined by equality (u 2 ) = T/Ga, where a is the lattice constant or 
the intermolecular distance. The average shear displacement is found from the relation- 
ship (u 2 ) = a 2 a 2 /G 2 . The total mean-square displacement represents the sum of these 
expressions provided that the thermal fluctuations and the stress are independent. Above 
implies that the transition of lubricant from solidlike to fluidlike state is induced both by 
heating and under influence of stress generated by solid surfaces at friction. This agrees 
with examination of solid state instability within the framework of shear and dynamic 
disorder-driven melting representation in absence of thermal fluctuations. Thus, the strain 
fluctuations, related to the stress ones, and the thermal fluctuations will be considered 
independently. 

In present paper the additive noises of the shear components of the elastic stress and strain, 
and the temperature are taken into account in a lubricant film. The solidlike lubricant is 
assumed to be amorphous (disordered). Therefore I study the glass transition represented 
in terms of a second-order phase transition. It is shown that increase of the strain noise 
intensity causes the lubricant film melting even at low temperatures of friction surfaces, 
at which the temperature noise plays a crucial role. The phase diagrams are calculated 
defining the domains of sliding, stick-slip, and dry friction in the planes temperature noise 
intensity — temperature of friction surfaces and noise intensity of shear elastic strain - 
temperature noise intensity. 



2 Langevin and Fokker-Planck equations 

Consider now the affect of additive noises of the elastic stress and strain shear components 
a, e, and the temperature T. With this aim, one should add to right-hand side of Eqs. (2) 
- (4) the stochastic terms I^ 2 ^, I^ 2 £, It 2 £, (here the noise intensities J CT , J e , It are 
measured in units of a 2 , e 2 r~ 2 , (T c k/1) 2 , correspondingly, and is the 5-correlated 
stochastic function) [19]. Then, within the adiabatic approximation r a ^> r e , tt, equations 
(3) and (4) are reduced to the time dependencies 



Here, deterministic components are reduced to obtained in Ref. [15], whereas fluctua- 
tional ones follow from the known property of variance additivity of independent Gaussian 
random quantities [19]. Thus, using the slaving principle inherent in synergetics [17,18] 
transforms noises of both strain e and temperature T, which are adiabatic initially, to 
multiplicative form. As a result, a combination of Eqs. (2), (5), and (6) leads to the 
Langevin equation 



e(t)=e + e^t), T(t) — T + T£(t); 



(5) 




f=(T e +2a 2 )d(a), T = J I T +I e a* d(a), d(a) = (1 + a 2 )' 1 



(6) 




(7) 



where the force / is related to the synergetic potential [15] 

V= 1 -(l- g y + g(l-^Hl + a 2 ) (8) 

and an expression for the effective noise intensity 

1(a) = I a + (j e + I T a 2 ) g 2 d 2 (a) (9) 

is obtained in accordance with above mentioned property of noise variance additivity. In 
order to avoid mistakes, one should notice that a direct insertion of Eqs. (5) and (6) into 
(2) results in the appearance of a stochastic addition 

il' 2 + {i l J 2 + 4 /2 -) gd{*)] m (io) 

whose squared amplitude is quite different from the effective noise intensity (9). Moreover, 
in contrast to the expressions (6), a direct use of the adiabatic approximation in Eqs. (3) 
and (4) reduces the fluctuational additions in Eqs. (5) to the forms: e = (I^ 2 + I^ 2 a)d(a), 
f = (4 /2 - Pj 2 a)d(a). The latter is obviously erroneous since the effective noise of the 



temperature T disappears entirely for the stress a = yI T /I £ . The reason for such a 
contradiction is that the Langevin equation does not permit the use of usual analysis 
methods (see [19]). 

To continue in the usual way, let us write the Fokker-Planck equation related to the 
Langevin Eq. (7): 

^) = ^{_ /W P ( ,, t)+ ^ [/W P ( ,, t)] }. 



At steady state, that is the single considered case, the probability distribution P(a, i) 
becomes a time-independent function P(a). Consequently, under the usual condition, 
that the expression in braces of the right-hand side of Eq. (11) is equal to zero, this leads 
to a stationary distribution 

P(a) = Z- 1 ex V {-U(a)}, (12) 
where Z is a normalization constant. The effective potential 



is determined by the synergetic potential V [Eq. (8)] and the noise intensity 1(a) [Eq. (9)] 
[20]. Combining these expressions, one can find the explicit form of U(a), which is too 
cumbersome to be reproduced here. The equation defining the locations of the maximums 
of the distribution function P(a) 



(1 - g)x 3 + g(2-T e )x 2 -2g 2 I T x+4g 2 (I T -I £ )=0, x = l + a 2 , 



(14) 



is much simpler. According to Eq. (14), maximums are insensitive to changes in the 
intensity of the noise I a of the stress a, but they are determined by the value T e of the 
friction surfaces temperature and the intensities I £ and I T of the noises of the strain e and 
the lubricant film temperature T, which acquire the multiplicative character in Eq. (9). 
Hence, for simplicity I a can be set equal to and Eqs. (8), (9), and (13) give the following 
expression for the effective potential: 



U(a) 



i -g)-g(2- T e )\ In |1 + a/i\ + (1 - g) 



a 



+ [g (2 - T e ) - i(l - g)} a} + In [g 2 d 2 (a) (l £ +a 2 I T ) 



i= — -l. 



(15) 



3 Phase diagrams 



According to Eq. (14), the effective potential (15) has a minimum at a = if the temper- 
ature T e does not exceed the critical level 



T c =l + g- 1 + 2g(I T -2I £ 



(16) 



whose value increases at increasing the characteristic value of shear viscosity i]q and the 
temperature noise intensity, but decreases with growth of the shear modulus G of lubricant 
and the strain noise. Here, the lubricant film does not melt. The solutions of Eq. (14) 



g(T e -2) 

1-9 



-3± 




(17) 



are obtained on the line defined by expression (16) after elimination of the root a 2 = 0. 
At T e <T c °=2(l+2g~ 1 ) /3 the roots a± are complex, starting from T e = T c ° they become 
zero, and at T e > T c ° one has real magnitudes a + = —o_ ^ that implies lubricant 
film melting. If equality (16) is fulfilled, the root o = corresponds to the minimum of 
the effective potential (15) at T e <T c °, whereas at T e > T c ° this root corresponds to the 
maximum, and the roots a± — to symmetrical minimums. 

Now, let us find another condition for the stability of the roots a± in the simple case 
I £ — 0. Setting the discriminant of Eq. (14) equal to zero, one gets the equations 



0, I 2 



27 
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8(1 -g) J 2g(l-g) 



0, 



the second of which gives 



2/ T = ^ 
2g 



'1-9 , 2-T e 



(2-T e ) 2 

8(1 -g) 



± 



27 (l-g 2-T e \ (2-T e ) 2 
2g{ g 3 J 8(l-<?) 



2 + 2 ( 2 - T -) 3 l V2 (19) 



This equation defines a bell-shaped curve T c (It), which intersects the horizontal axis at 
the point 



9(3-4fl)+8fl 2 
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(20) 



and vertical axis at the point T e = 2. It has a maximum T e = 2g 1 at 
r 2(1-^) 



(21) 



It is easy to see that line (16) touches the curve (19) at the tricritical point 



(22) 



Thus, this point addresses to the appearance of real roots a± ^ (17) of Eq. (14) that 
means lubricant film melting. 

Let us now consider the more general case of two multiplicative noises I £ , It ^ 0. The 
condition of extremum of the effective potential (15) splits into two equations, one of 
which is simply a = 0, and the other one is given by Eq. (14). As mentioned above, the 
analysis of the latter indicates that the line of existence of the zero solution is defined by 
expression (16). The tricritical point T has the coordinates 



T e =^(l+2( ? - 1 -2^/ £ ), / T =l(^- 1 -l+8^/ £ ). 
3 6g 



(23) 



The phase diagrams for the fixed intensities I £ are shown in Fig. 1. Physically, one should 
take into consideration that lines 1 and 2 define the thresholds of stability loss of the 
system. Above straight line 1 the system manifests a stable sliding friction (SF) inherent 
in the liquidlike phase of lubricant, below curve 2 the dry friction (DF) occurs that is 
characteristic for the solidlike state of lubricant film. Between these lines the region of 
the stick-slip friction (SS) mode is realized, i.e., mode that is characterized by periodic 
transitions between two dynamic states during steady-state sliding. It is relevant to an 
intermittent regime of lubricant melting, where a mixture of both solidlike and liquidlike 
phases exists. For I £ < (1 + g~ l )/Ag the situation is generally the same as in the simple 
case I £ — (see Fig. la). At I £ > (1 + g _1 )/4g the sliding friction is possible even for 
small values of temperature T e of friction surfaces and noise intensities It of the lubricant 
film temperature (Fig. lb). According to (23), the tricritical point lies on the Jy-axis at 
I £ = (1 + 2g~ 1 )/2g ) and if the noise intensity I £ is larger than the critical value I £ =2g~ 2 , 
the stable dry friction domain disappears (see Fig. lc). It is worth noting that this domain 
decreases with increase of the shear modulus value G and decrease of the characteristic 
value of shear viscosity i]q. 




Fig. 1. Phase diagrams at 5=0.5 and fixed values I £ : (a) I £ =l.2; (b) I £ =5; (c) I £ =8. Lines 1 and 
2 define the boundary of stability domains of sliding (SF), dry (DF), and stick-slip (SS) friction 
(T is the tricritical point, C is the critical point). 

The consideration of the additive noises of a, e, and T shows that the stochasticity 
influence is non-essential for the shear component of elastic stress tensor and it is crucial 
for both the corresponding component of strain and the lubricant film temperature. The 
boundary of the domain of sliding friction is fixed by the equality for the noise intensities 



I T = 2I e -(l+g)/2g 2 , 



(24) 



following from Eq. (14) at the conditions x—1 (cr—0), and T e =0. According to Eq. (24), 
in absence of the temperature noise the lubricant melting occurs if the noise intensity of 
the shear strain component exceeds the value 

I £0 = (l+g)/4g 2 , (25) 



corresponding to the point O in Fig. 2. The increase of both the shear strain and the 
temperature noises causes the lubricant melting if their intensities are bounded by con- 
dition (24). The domain of the stick-slip friction appears with further increase of these 
intensities above magnitudes 

/ - {2+9) I - 3+9 (26) 



at the tricritical point T in Fig. 2. Such an intermittent behavior is realized within the 
region located above straight line (24) and outside the curve that is determined by 
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If the noise intensity of the shear strain exceeds the value I e i defined by (27) with the tem- 
perature noise 1^2=2(3— g)/g 2 (the point C in Fig. 2), the dry friction region disappears 
at all. The curve (27) intersects the vertical axis at the point 
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(28) 



above that the dry friction does not take place. The corresponding phase diagram depicted 
in Fig. 2 has a very non-trivial form (especially, within the domain I e i<I e <I e 2)- 



4 Summary 



The above consideration of the thermal and elastic fields noise influence on the solid-liquid 
transition of ultrathin lubricant film permits to define the domains of dry, sliding, and 
stick-slip friction modes in the phase diagram. So that, an evidence of the phase diagram 
complication is obtained due to studied fluctuations. Depending on the initial conditions 
the growth of lubricant film's temperature noise can decrease or increase friction, but the 
growth of elastic shear strain noise increases the sliding friction region only. It is shown 
that dry friction domain is bounded by relatively small values of the confining walls tem- 
perature and the noise intensities of lubricant strain and temperature (see Figs. 1 and 2). 
Thus, used here approach predicts the possibility for controlling of frictional behavior. 




Fig. 2. Phase diagram for the system with T e =0, 5=0.5, and I £ , It / 0. 

Above the concept of dynamical shear melting of the ultrathin lubricant film has been 
used [10,11]. In accordance with it the stick-slip friction can be described and such melt- 
ing is represented as a result of action of elastic field of shear stress component caused 
by the heating of friction surfaces above the critical value. The essential limitation of 
this approach is the fact that stick-slip motion I studied was independent of the way in 
which the system was driven, i.e., elasticity and mass of the confining walls, although for 
such friction mode the mentioned dependence is crucial. It is worth noting that here the 
temperature of the confining walls T e plays a role of the parameter of external influence. 



Besides, the friction force decreases with increasing velocity at the contact V = Ide/dt 
because the latter leads to the growth of the shear stress a according to the Maxwell-type 
stress - strain e relation: da/dt— —a/r^+Gde/dt. 
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